Propagating, evanescent, and localized states in carbon nanotube-graphene junctions 
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We study the electronic structure ol the junctions between a single graphene layer and carbon 
nanotubes, using a tight-binding model and the continuum theory based on Dirac fermion fields. The 
latter provides a unified description of different lattice structures with curvature, which is always 
localized at six heptagonal carbon rings around each junction. When these are evenly spaced, we 
find that it is possible to curve the planar lattice into armchair (6n, 6n) as well as zig-zag (6n, 0) 
nanotubes. We show that the junctions fall into two different classes, regarding the low-energy 
electronic behavior. One of them, constituted by the junctions made of the armchair nanotubes 
^ , and the zig-zag (6n, 0) geometries when n is a multiple of 3, is characterized by the presence of 

two quasi-bound states at the Fermi level, which are absent for the rest of the zig-zag nanotubes. 
These states, localized at the junction, are shown to arise from the effective gauge flux induced 
CNj , by the heptagonal carbon rings, which has a direct reflection in the local density of states around 

the junction. Furthermore, we also analyze the band structure of the arrays of junctions, finding 
out that they can also be classified into two different groups according to the low-energy behavior. 
In this regard, the arrays made of armchair and (6n, 0) nanotubes with n equal to a multiple of 3 
are characterized by the presence of a series of flat bands, whose number grows with the length of 
the nanotubes. We show that such flat bands have their origin in the formation of states confined 
to the nanotubes, with little overlap in the region between the junctions. This is explained in 
. the continuum theory from the possibility of forming standing waves in the mentioned nanotube 

geometries, as a superposition of modes with opposite momenta and the same quantum numbers 
under the Cev symmetry of the junction. 
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I. INTRODUCTION 



^ , The recent isolation of graphene layers a single atom thickfT], 0, [1] has lead to a great deal of activity, because 
Q ' of their novel electronic properties and potential applications. The lattice structure of graphene is determined by 
^ , the sp^ coordination between neighboring carbon atoms. Each carbon atom has three nearest neighbors, leading to 
planar honeycomb lattice. With small modifications, the same structure describes other carbon allotropes, such as 
the fuUerenes and the carbon nanotubes. 
^ . In this paper we study the simplest systems which combine two of these allotropes: the junctions between a single 
f — ' graphene layer and carbon nanotubes. A regular array of closely spaced armchair nanotubes attached to a graphene 
lO layer has already been studied 3, and related systems are being considered for their potential applications 5]. We 
undertake here the investigation of junctions made of nanotubes with different chiralities, which have in common a 
^""^ transition from the planar to the tubular geometry mediated by the presence of six heptagonal carbon rings. These 
T-H ' induce the negative curvature needed to bend the honeycomb carbon lattice at the junction, playing a kind of dual 

I role to that of the pentagonal carbon rings in the fullerene cages @]. 
0^ We analyze first the electronic properties of a single junction between a carbon nanotube and a graphene layer. 

We discuss the possible structures of this type, concentrating on geometries where the heptagonal rings are evenly 
spaced around the junction. The nanotubes can be then either armchair (6n, 6n) or zig-zag with (6n, 0) geometry 
(that is, with 6n hexagonal rings around the tube). We calculate their electronic structure, using the tight-binding 
^ , model based on the tt orbitals of the carbon atoms widely applied to carbon allotropes with sp^ coordination. 

' Paying attention to the local density of states, we find that the junctions fall into two different classes, depending 
on the behavior in the low-energy regime. One of the classes, comprising the junctions made of armchair and (6rt, 0) 
nanotubes when n is a multiple of 3, is characterized by the presence of a peak in the density of states close to the 
Fermi level. The peak is absent in the other class, formed by the junctions made with the rest of zig-zag geometries. 
In general, the density of states tends to be depleted in the junction at low energies, with peaks above and below the 
Fermi level marking the threshold for the propagation of new states across the junction. 

We present next a continuum description, based on the formulation of Dirac fermion fields in the curved geometry, 
which allows us to characterize the general properties of the junction, and which is consistent with the previous discrete 
analysis. Thus, we see that the peak at the Fermi level in the local density of states is in general a reflection of the 
existence of quasi-bound states (zero modes) for the Dirac equation in the curved space of the junction. It is known 
that the topological defects of the honeycomb lattice (pentagonal and heptagonal rings) induce an effective gauge 
field in the space of the two Dirac points of the planar graphene lattice [3]. It turns out that the effective magnetic 
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flux is enough to localize two states at the junctions made of armchair or (6n, 0) nanotubes when n is a multiple of 3. 
At low energies, however, the generic behavior is given by evanescent states, which arise from the matching of modes 
with nonvanishing angular momentum and have exponential decay in the nanotube. 

We finally apply our computational framework to the analysis of the band structure of the arrays of nanotube- 
graphcnc junctions. Considering the behavior of the low-energy bands close to the Fermi level, we find that the arrays 
also fall into two different classes. The arrays made of armchair nanotubes or (6n, 0) nanotubes with n equal to a 
multiple of 3 tend to have a series of flat bands close to the Fermi level, while the arrays made with the rest of zig-zag 
nanotubes have all the bands dispersing at low energies. Such a different behavior has its origin in the existence of 
states confined in the nanotube side of the junction. We find that this feature can also be explained in the context of 
the continuum model. The armchair and the (6n, 0) geometries with n equal to a multiple of 3 allow for the formation 
of standing waves between the junction and the other end of the tube. This is the mechanism responsible for the 
confinement of the states in the nanotubes and the consequent development of the flat bands, whose number grows 
at low energies with the length of the nanotube, in agreement with the predictions of the continuum theory. 

II. TIGHT-BINDING APPROACH TO CARBON NANOTUBE-GRAPHENE STRUCTURES 

A. Lattice structure 

Our first aim is to analyze the density of states of a semi-infinite nanotube attached to a graphene layer in the 
tight-binding approximation. The possible setups that we will consider, keeping the threefold coordination of the 
carbon atoms, are sketched in Fig. [1] The structures can be wrapped by the graphene hexagonal lattice, with the 
exception of the six points where the sides of the hexagonal prism (which describes the nanotube) intersect the plane. 
The threefold coordination of the carbon atoms requires the existence of sevenfold rings at those positions. 




(a) (b) (c) 

FIG. 1: (a) Sketcli of carbon nanotube attached to a graphene plane. The building blocks of the structure are triangles which 
include many carbon atoms. The orientation of the bonds of the honeycomb lattice in the basic triangle may give rise to 
armchair (b) or zig-zag nanotubes (c). The threefold coordination of the carbon atoms induces the existence of six heptagonal 
rings at the ends of the sides of the hexagonal prism contacting the plane in (a). 

We describe the electronic states in the structures shown in Fig. [1] by means of a nearest-neighbor tight-binding 
model. In general the relaxation of elastic energy will modify the bond lengths at the junction, depending on the 
nanotube radius. We will assume that this relaxation does not change significantly the electronic behavior. In this 
respect, a tight-binding model based on the tt carbon orbitals is well-suited for the purpose of discerning the extended 
or localized character of the different electronic states. Our main achievement will be to assign the different features 
in the local density of states to the behavior of the electronic states near the nanotube-graphene junctions. To this 
aim, we have actually checked that slight modulations of the transfer integral t near the junctions do not produce 
significant changes in the results shown in what follows. 

B. Electronic densities of states 

We concentrate on the analysis of geometries where the six heptagonal carbon rings are evenly spaced around the 
junction as in Fig. [1] This constrains the possible chiralities of the nanotubes, that can be then either armchair 
(6rt, 6n) or zig-zag (6n, 0), with the number n running over all the integers. Nanotubes in which the carbon sheet is 
wrapped with helicity can be also attached at the expense of introducing an irregular distribution of the heptagonal 
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FIG. 2: Sequence of local densities of states for a circular ring of atoms at the end of a (54, 0) nanotube close to the junction, 
for the different sectors corresponding to eigenvalue q under n/3 rotation equal to 1 (a), e**'^''"^ (b), e*^"^/"^ (c), and —1 (d). 
Energy is measured in units of the transfer integral t. 



rings. Anyhow, we expect that the rules explaining the different features in the density of states are universal enough 
to hold even in these more general cases. 

We have obtained the spectra of different types of hybrid structures by diagonalization of the tight-binding hamil- 
tonian for very large lattices, with up to « 50, 000 carbon atoms in the graphene part and « 40, 000 in the nanotube 
side. Given that the whole geometry has Cq^ symmetry, we have classified the energy eigenstates into six groups 
according to the eigenvalue q under a rotation of 7r/3. The nature of each electronic state is given in general by its 
behavior at the nanotube-graphene junction. For this reason, we have characterized the hybrid structures in terms of 
the local density of states averaged over a circular ring of atoms at the end of the nanotube close to the junction. 

Our computations have covered a number of structures including armchair and zig-zag nanotubes with different 
radii. After inspection of all the spectra, it becomes clear that there are several generic features in the density of 
states. We have represented in Fig. [2] the behavior near the junction between a graphene layer and a (54,0) zig-zag 
nanotube. We observe that, apart from the peak close to zero energy in the sectors corresponding to q = e^"^^^^, 
for q =/= 1 there is always a depletion in the density of states at low energies, delimited by two abrupt upturns. It is 
remarkable that the pattern in the sectors corresponding to q = q'^'^™!'^ reproduces the same observed for q = e^™/^, 
but with a scale that is approximately twice larger. The density of states in the sector with q = —1 displays in turn 
a wider depletion, with the position of the peaks scaled by an approximate factor of 3 with respect to those in the 
q = gi^'^/s sectors. 

In the above behavior of the density of states, the appearance of the peak close to zero in the sectors with q — 6=*=*^/^ 
is the only feature not generic for all kinds of nanotubes. In the junctions made with zig-zag nanotubes, the peak 
actually appears for nanotube geometries of the type (6n, 0) when n is a multiple of 3. In this series of hybrid structures, 
the patterns in the density of states for each value of q are quite similar, with the position of the corresponding peaks 
scaled in proportion to the radius of the nanotube. On the other hand, the rest of junctions, for which n is not a 
multiple of 3, display a different behavior. We have represented in Fig. [3] the local density of states averaged over a 
ring of atoms at the end of a (48, 0) nanotube close to the junction. It can be observed the depletion of the density 
of states at low energies in all but one of the g-sectors, and the absence of a peak at zero energy in any of the sectors. 

The local density of states of the (48, 0) nanotube is dominated at the junction by contributions from states with 
q = e±2i7r/3^ g^-^^ ^Yas has to do with the fact that the lowest-energy subbands of a (6n, 0) zig-zag nanotube have a 
nonvanishing angular momentum equal to ±4n for the motion around the tubule. This corresponds to a quantum 
number q = e^"^'"/^ in the case of the (48, 0) nanotube, while the low-energy states have g = 1 in the (54, 0) nanotube. 
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FIG. 3; Similar sequence as in Fig. [2] for a (48, 0) nanotube close to the junction, for q = 1 (a), e^^""'^ (b), e^^'''''^ (c), and -1 
(d). 

The present picture becomes then consistent with the fact that states in higher subbands may propagate across the 
junction only above (or below) some threshold energy. This feature will be established more precisely in the continuum 
approach derived below in terms of the Dirac equation. 

At this point, the hybrid structures can be classified into two different groups, depending on whether there is a peak 
or not close to zero energy in the local density of states around the nanotube-graphene junction. The peak comes 
actually from the contribution of a doubly degenerated level with states having q — e^"/^, and whose probability 
distribution decays exponentially in the nanotube. The character of these states will be established in the next section, 
after developing the continuum limit in terms of Dirac fermion fields. 

The evidence for the two different classes of hybrid structures is reinforced by the fact that the junctions made with 
armchair nanotubes behave in a quite similar way to that shown by the (6n, 0) nanotubes when n is a multiple of 3. 
The local density of states averaged over a circular ring of atoms around the junction between a (12, 12) nanotube 
and a graphene layer has been represented in Fig. |4l We observe the presence of the peak close to zero energy in the 
sectors with q = e^™!'^ . There is a clear depletion in the local density of states at low energies except in the sector 
with q — \^ which is consistent with the fact that the lowest-energy subbands in the armchair nanotube correspond 
to zero angular momentum around the tubule. 

The existence of the two different classes of nanotube-graphene junctions is illustrated in Fig. [5l which shows the 
results for the local density of states around the junction (after summing over the sectors with different values of q) 
for the different types of nanotube considered above. It is remarkable the similarity between the density of states for 
the (18,0) and (12,12) nanotube geometries, which have a very close value of the radius. This suggests that there 
must be a universal way of understanding the low-energy electronic properties of the two different classes of junctions, 
independent of the details of the lattice building the junction within each class. 



III. CONTINUUM APPROACH 



A. Heptagonal rings, effective flux, and matching at the junction 



We observe that the density of states summed over all values of q tends to have an approximate linear behavior 
away from the very low energy regime. This motivates the analysis of the junction in terms of the Dirac equation in 
the hybrid geometry, as the Dirac fermions provide an appropriate description of the electronic properties in the ID 
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FIG. 4: Similar sequence as in Fig. [2] for a (12, 12) nanotube close to the junction, for g = 1 (a), e=^'''/^ (b), e^^'""^^ (c), and 
-1 (d). 




FIG. 5: Sequence of local densities of states for a ring of atoms at the end of the nanotube close to the junction, for the different 
geometries (54,0) (a), (18,0) (b), (48,0) (c), and (12, 12) (d). Energy is measured in units of the transfer integral t. 
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carbon nanotube as well as in the 2D graphene layer. 

We will then assume that the radius i?o of the nanotube is much larger than the graphene lattice constant, in order 
to obtain the continuum limit of the tight-binding model. Within this approximation, the analysis of the electronic 
structure is reduced to the study of the Dirac equation in a space with an abrupt change from a planar to a cylindrical 
structure. The transition from a geometry to the other takes place due to the presence of the six heptagonal rings 
at the junction. These defects are the source of negative curvature, playing a role opposite to that of the pentagonal 
rings in a fuUerene cage ^6]. 

The heptagons, as well as the pentagons, also induce frustration in the honeycomb lattice, and lead to the exchange 
of the two Dirac valleys of the planar geometry Q. This latter effect has to be accounted for by means of an effective 
non-abelian gauge field operating in the space of the two independent Dirac points of graphene. It can be shown 
that the effective flux associated to an individual heptagonal ring is equal to 7r/2 [7] (in units such that h = 1). The 
flux provided by the six heptagons at the junction can reach therefore a maximum of 3/2 times the flux quantum. 
In general, however, the count of the total flux may not follow an additive rule, so that it can be lower than the 
maximum value, depending on the relative position of the heptagonal ringsQ. We will see that this is actually the 
origin of the two different classes of junctions. 

Note that the angular momentum around the axis of the nanotube is conserved and quantized in integer units, 
in a continuum description of the geometry analyzed here. On the other hand, the angular momentum in the 
plane is quantized and shifted by 1/2 plus the number of flux quanta of the effective gauge field induced by the 
heptagonal rings. The existence of topological defects which induce an effective flux at the junction allows us to 
match wavefunctions with different angular momenta at either side of the junction, provided that the effective flux 
corresponds to a half-integer number of quanta. 



B. Continuum wavefunctions 



By looking at the effect of pairs of heptagonal rings, the non-abelian gauge field operating in the space of the two 
Dirac points becomes anyhow proportional to a sigma matrix T3. We may consider then that the effect of the six 
heptagonal rings at the junction is described by an effective abelian field, standing for either of the eigenvalues of T3 
0. We end up therefore with two different Dirac equations, with effective magnetic fluxes of opposite sign. We will 
denote the Dirac spinors satisfying the two Dirac equations as 5*+ and respectively. Given that the region of 
the space away from the junction has no curvature, we can use radial coordinates r, 9 to write the set of two Dirac 
equations in the plane for r > Rq: 

e^^ir^e) (1) 

e^iirJ) (2) 

vp the Fermi velocity and the components '^a and 'i's denote the respective amplitudes of the electron in the two 
sublattices of the graphene lattice. The parameter g corresponds to the quanta of effective flux felt by the electrons 
when making a complete tour around the junction. 

In the nanotube side, we use cylindrical coordinates z, 9, with z < 0. The Dirac equation for the nanotube is 

e^^{z,9) (3) 

e^^{z,9) (4) 

We note that (HJ-lHl) as well as ©-([I]) are expressions of the Dirac equation in flat space. This is consistent with 
the fact that, in the continuum limit, the curvature is localized at the circle connecting graphene and the nanotube. 
As an alternative to the coordinate z, we could make for instance the change of variables r — Rq exp(z/i?o), allowing 
us to map the nanotube into the region of the plane with r < Rq. Using a common radial coordinate r > to 
describe both graphene and the nanotube, the metric of the space turns out to be multiplied by the conformal factor 
r2(r) = 9{r — Rq) + (i?o/?')^6'(i?o — r). The first derivative of the metric becomes then discontinuous and the curvature 
scalar, computed in terms of second derivatives of the metric, is i? = —2rS{r — Rq). The integral of this expression 
corresponds actually to the total curvature provided by the heptagonal rings. This makes clear that the effects of the 
curvature are implicit in the operation of matching the solutions of ([ij-lll]) and ©-([l]) at the circle r — Rq. 



Rq 
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Rq 
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The solutions of the Eqs. JlJ-llll) are of the form 

where ci and C2 are constants, and Jn{x) and Yn{x) are Bessel functions. The energy is e = ivpk. 

On the other hand, the resolution of Eqs. ©-([l]) shows that there are propagating and evanescent waves in the 
nanotube, which can be written as: 



1 \ ... , / 1 



vl/±(z, e) \ _} "=11^ -sgn(e) e^^C^) ]^ + sgn(£) e'^^^^') )^ ^ ^ 



-sgn(e) e^'^^-''") / ' ' " «o 

where c'j^, C2 and c are constants, the energy e is given by: 



(6) 



and the phase factor in Eq. ^ is: 

gi0(fc) ^ (8) 



%n 



We note that the evanescent states with longitudinal decay e'^^ arise for nonvanishing angular momentum n. Then 
there is an energy threshold VF\'n\IRQ for the appearance of propagating states in the nanotube. This is perfectly 
consistent with the behavior of the local density of states obtained for the different values of q in the tight-binding 
approach. The depletion found in different g-sectors for the local density of states at the end of the nanotube (close to 
the junction) corresponds actually to the range of evanescent states given by (l7|). As already mentioned, the position 
of the peaks delimiting the depletion in the tight-binding approach scales in proportion to the value of the angular 
momentum, which corresponds in the lattice to the different values of q. Moreover, we have also seen that such 
a position is inversely proportional to the nanotube radius _Ro, with values in the plots that can be approximately 
matched with the estimate vftl/Rq (after using the expression of the Fermi velocity vp — 3to/2, in terms of the 
transfer integral t and the C-C distance a). We find therefore that the generic features found for the local density of 
states in the tight-binding approach are well captured by the continuum limit based on the Dirac equation. 



C. Zero-energy states 

The presence of a peak in the local density of states close to zero energy (in the sectors q — e^^^^/^) jg the only 
feature not generic for all types of nanotubes, and that can be also explained within our continuum approach. The 
rotation caused by each heptagonal ring in the space of the two Dirac points corresponds to an effective magnetic 
flux of 7r/2 0, but the way this flux is combined in the case of pairs of heptagons depends on their relative position. 
This has been studied in the case of pentagon pairs in Ref. arriving at a conclusion that can be readily generalized 
to the case of heptagonal defects. The result is that, when the distance between the heptagons is given by a vector 
(iV, M) (using the same notation to classify the nanotubes) such that — Af is not a multiple of 3, the effective flux 
of a pair of heptagons does not add to tt, but to the lower amount tt/S. The number N — M for the distance between 
heptagons is a multiple of 3 only in the case of junctions with armchair nanotubes, or with (6n, 0) nanotubes when n 
is a multiple of 3. In these instances, the total flux felt around the junction is equal to the sum of the fluxes provided 
by the individual heptagons, giving a value of g = 3/2. In the rest of the cases, the total flux corresponds instead to 
5 = 1/2. 

The number g of flux quanta has a direct correspondence with the number of zero modes of the Dirac equation. 
Their existence rests on the possibility of having localized states at the junction, with suitable decay in both the 
graphene part and the nanotube side. If we take for instance the maximum effective flux and g — 3/2, we have an 
equation for a zero-energy eigenstate in the region r > Rq 



\ r r J 



(9) 
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For a wavefunction with angular momentum n, we obtain the behavior 

^A{r,e) r^r'^-^e'"^ (10) 

This gives rise to modes decaying from the junction for values n < I . The wavefunction has to be matched at r = Rq 
with the appropriate dependence in the nanotube, that is 

^'a(z,6') - e^^e"* (11) 

Recalling that z < 0, we see that only the value n = 1 provides a localized state at the junction. On the graphene 
side, the state is not strictly normalizable, in a similar way to other half-bound states induced by defects[lll, [T3 |. 
Anyhow, such a localized state has a reflection in the peak observed close to zero energy in the q — e^"/^ sectors of 
the tight-binding density of states. By inverting the direction of the flux and taking g — —3/2, it can be seen that 
the solutions have then a nonvanishing component ^'_B(r, 0) similar to ()10|1 . but with angular momentum —n instead 
of n. Another localized state is found therefore with opposite chirality and n = —1. 

In the case of the junctions with (6n, 0) nanotubes such that n is not a multiple of 3, the flux corresponding to 
g — 1/2 is not enough to localize states at the junction. It can be seen that there are no zero-energy solutions of the 
Dirac equation decaying simultaneously in the graphene plane and in the nanotube. This explains why in this type 
of junctions there is no low-energy peak within the depleted region of the local density of states. We complete in this 
way the correspondence between the tight-binding approach and the continuum limit based on the Dirac equation, 
accounting for the main electronic features and unveiling also the origin of the different low-energy behavior in the 
two classes of junctions. 



D. Transmission into the nanotube 



Eqs. ^ and ^ allow us to analyze the scattering of a wave in the plane off the nanotube. The coefRcients ci, C2, c'l 
and c'2 define the transmission and reflection by the nanotube. For a wave coming from the plane, the coefficient c'2 
is zero, and c'^ is proportional to the transmission coefficient. Using the theory of scattering of two dimensional Dirac 
electrons by an impuritvfol. [isl. [3. [Tsl. [iGj . the transmission coefficient is given by: 



7rfc_Ro iJn+iikRo)Yn{kRQ) - iJn{kRo)Yn+i{kRo) 
2 ir„(fci?o) -r„+i(fci?o)e'*W 



At high energies, e S> vp/Rq or, alternatively, fci?o ^ 00, we can use the asymptotic expansion: 



(12) 



and hmfc_Ro 



lim Jn{kRo) , , 

kRa^oo y TTfcilo 

lim y„(fci?o) « 
1. From these expansions, we obtain: 



, TiTT vr 

COS ( kRo - - ^ 

, nil TT 
sm I fci?o 2~ ^ 4 



hm r„ 

/c_Ro — 'OO 



zsin fci?o-^-f 



cos 



i[kRQ-^~l) 



(13) 



(14) 



so that limfc/jjj^oo \Tn\'^ = 1- This estimate is valid for angular momenta n such that n <^ fci?o- We can also obtain 
the reflection coefficient in this limit, Rm which is also independent of n for n <C fc-Ro- The angular dependence of 
the scattering cross section a{0) is: 



sin{kRo9) 
sm{e/2) 



(15) 



and the total cross section, a — J a{9)d9, is proportional to Rq. The total flux of particles propagating inside the 
nanotube, normalized to the total incoming flux, is also proportional to /Jo- 
in the low energy limit, fci?o ^ 1, we obtain: 



lim |T„|^ 

kRo^O 



kRo < n 



(16) 



In this limit, most of the electrons reaching the junction are scattered back into the plane. 
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E. Local density of states 



We can make use of the continuum equations to analyze systems of very large sizes. We calculate the electronic 
Green's functions numerically. The Dirac equation in the plane, in radial coordinates, can be discretized[l7|. Each 
radial equation can be approximated by a nearest neighbor tight binding model with two inequivalent hoppings: 



2n+l 
2n+ 1 



2n+ 1 



2n+ 1 



4i 



b, 



eh. 



eai 



(17) 



where and hi give the values of ^'^(r) and at position r = i x a, a being a length scale which defines the 

discretization. Eqs. (|17p also include an energy scale, t, which plays the role of an upper cutoff. The Fermi velocity 
is equal to t x a. The Dirac equation is obtained for e <^t. We can write Eqs. (|17p in a more compact form using a 
single index, such that C2j-i — aj and C2j — bj, so that: 



tiCi-i + ii+iCi+i — eci 



where, using as new length scale a/2, we have: 



t. = t 



i-i-iy 



2n+l 



4(z + 7Vo)_ 

where we start at position A^o- The diagonal Green's function at site i can be written as: 

1 



and: 



with boundary conditions at i — and i — N: 



Giiie) = 



e - UTl 



1 



1 

e - SAr(e) 



and: 



(18) 



(19) 



(20) 



(21) 



(22) 



So(e) 
Sjv(e) 



<^ 1 



2No 



it 



and: 



t 



2 



v/?2 - 4t2 



e = e-- 1 



{n-lf 



(23) 



(24) 



where the radius of the nanotube is Rq = Nq x a/2. The boundary conditions in Eq. (j23p describe a semi- infinite 
nanotube attached at position Nq, and also approximate the boundary of a plane at position N. The Green's function 
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FIG. 6: Density of states at the bulk of the nanotube (blue, hexagons), and at the graphene plane (red, diamonds), obtained 
using the numerical method discussed in the text. 
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FIG. 7: Density of states at a graphene plane attached to a nanotube. The radius of the nanotube is A'o = 100. The density 
of states is shown at position n = 10 from the junction (red, diamonds), and n = 100 (blue, hexagons). 

deep inside the nanotube can be calculated analytically: 

The density of states in the plane has been calculated numerically, with N = 1600 and summing angular momenta 
from n — —100 to n = 100. The calculation is equivalent to analyzing a cluster with (2 x 100 + 1) x 1600 — 321600 
sites. Results for the Green's function in the nanotube, and at the position i = 100, are shown in Fig. [6] The energy 
scale is set by i = 1. 

The Green's function for the system built up by the nanotube and the graphene sheet is shown in Fig. [71 The 
parameters are the same as the ones used for the calculation shown in Fig. [51 and the radius of the nanotube is 
No = 100. 

The density of states at distances from the nanotube n > iVo are similar to those in the unperturbed sheet. Near 
the juncture with the nanotube, there is a depletion of states at low energies, compensated by the existence of a 
localized state at e = 0. 



F. Gauge fields due to elastic strains 

We have not analyzed so far the effect on the electronic structure of strains which may be induced near the junction. 
These strains deform the bonds, and induce an additional, intravalley gauge field acting on the electrons [3. ITsl [T9I. |20| . 
The large in-plane stiffness of graphene implies that the bonds will tend to their equilibrium lengths throughout the 
system. 
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The bending at the junction will be localized within a length scale I \/ k/A, where k ~ 1 eV is the bending rigidity 
of graphene, and A ~ 10 eV is an average of the Lame coefficients of graphene. This length is comparable to the 
lattice spacing. 

The mismatch between the diameter of the nanotube and the lattice constant of the graphene layer induces additional 
strains, with a long range decay into the bulk of the graphene plane and the nanotube, which can be calculated using 
the continuum theory of elasticity 2l|. We expect, however, the interatomic distance in graphene to be very close to 
the distance between carbon atoms along the radial direction of the nanotube, so that the strains induced by this 
effect will be small. The strains will decay as or as function of the distance to the junction. 

Using dimensional analysis, the strain near the junction is of order ARq/Rq, where Ai?o is the change in the 
equilibrium radius of the nanotube induced by the plane or, alternatively, of order Aa/a, where a is the interatomic 
distance. We expect the value of Aa/a to be similar on the plane side of the junction. The associated gauge field is 
A - /3Aa/a2, where /3 ~ i91og(i)/aiog(a) - 2 — 3 gives the change of the tight-binding hopping t with a. Thus, we 
expect that the elastic strains will induce changes on the electronic structure on energy scales of order vpA near the 
junction. 



IV. LOW-ENERGY BANDS IN ARRAYS OF NANOTUBE-GRAPHENE JUNCTIONS 

Our computational framework allows us also to address the electronic properties of arrays of nanotube-graphene 
junctions. We consider the case in which the unit cell of the array has a hexagonal shape in the base, of the type shown 
in Fig. [TJ The periodic arrangement of junctions is formed then by translating the unit cell by two independent vectors 
of the triangular array, in such a way that the 2D base is fully covered with the hexagonal patches. The Brillouin zone 
of the superlattice is a hexagon, and the main electronic properties are encoded in the form of the bands from the 
center to the M and K points at the boundary of the zone. As long as the states in momentum space have well-defined 
transformation properties under translations by the lattice vectors of the triangular array, we can obtain the bands 
of the array of junctions by solving a tight-binding model in the unit cell, with appropriate momentum-dependent 
boundary conditions between opposite sides of its hexagonal base. 

The band structure of the array of junctions depends on the geometry of the nanotubes, as well as on their length 
and the distance between them. For simplicity, we are going to consider arrays where all the nanotubes have the 
same chirality. Then, it can be checked that the arrays fall into two main classes, regarding the behavior of the bands 
close to the Fermi level. The distinctive feature of one class with respect to the other is the presence of flat bands 
in the low-energy part of the spectrum. The arrays of junctions made of armchair nanotubes, for instance, always 
have a number of these fiat bands, as illustrated by the representative in Fig. EJa). The appearance of fiat bands 
in a particular array of junctions was noticed in Ref. d. We have found that the flat bands are actually generic in 
arrays made of armchair nanotubes, which display a series of them as one moves from the Fermi level to higher (or 
lower) energies. The spacing in energy between the fiat bands becomes inversely proportional to the length of the 
nanotubes. The bands dispersing at low energies in Fig. [SJa) are not affected however by variations of that variable, 
while they move instead closer to the Fermi level as the distance between the nanotubes in the array is increased. 

On the other hand, the presence of fiat bands at low energies is not generic in arrays made of zig-zag nanotubes. 
In general, we may expect a number of bands dispersing above and below the Fermi level, as shown in Figs. [H^b) and 
(d), which represent the low-energy bands of arrays made respectively of (12,0) and (24,0) nanotubes. In the case 
of the zig-zag nanotubes, flat bands only appear close to the Fermi level when the junctions are formed with (6n, 0) 
geometries such that n is a multiple of 3. This distinctive behavior can be appreciated in Fig. [5{c), which displays 
the low-energy bands in the case of an array made of (18, 0) nanotubes. The shape of the bands resembles there the 
typical appearance of the spectra of arrays made of armchair nanotubes, as shown in Fig. [SJ^a). 

The mentioned fiat bands have their origin in the existence of localized states in the arrays of junctions. We have 
checked that the junctions made of armchair nanotubes and (6n, 0) nanotubes with n equal to a multiple of 3 have in 
common the formation of electron states confined mostly in the nanotube side. These are the states responsible for the 
development of the flat bands shown in Fig. [51 as the wave functions with most of their weight in the nanotubes show 
little overlap in the graphene part of the lattice. This also explains in a natural way the proliferation of fiat bands at 
low energies as the nanotube length is increased, by thinking of the confined modes as standing waves in the tube. On 
the other hand, the states that are preferentially localized at the junctions (corresponding to the quasi-bound states 
of the Dirac equation) may be identified here as the pairs of branches degenerated at the F point. These states may 
have in general a significant overlap between nearest junctions, which is reflected in the appreciable dispersion of the 
corresponding bands. 

We can reach in the continuum limit a qualitative understanding of the similar behavior of the arrays made of 
armchair and (6n, 0) nanotubes when n is a multiple of 3, by noticing that these are the only geometries that support 
low-energy standing waves between the junction and the other end of the tube. This requires the superposition of 
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FIG. 8: Low-energy bands of arrays of junctions with nanotube geometries (12, 12) (a), (12, 0) (b), (18, 0) (c), and (24, 0) (d). 
In all the cases the unit cell of the array is of the type shown in Fig. [T] with a side of the hexagon in the basal plane equivalent 
to 10 carbon rings of the graphene sheet, and a nanotube height equivalent to 10 unit cells of the nanotube for (a), (b), (c), 
and 20 nanotube unit cells for (d). 

two modes with opposite momenta along the tube, which is possible at low energies in the armchair nanotubes as the 
modes at opposite Dirac points have then vanishing angular momentum. In general, this is not the case for the (6n, 0) 
geometries, since in the zig-zag nanotubes the Dirac points correspond to large momenta in the transverse direction. 
Yet the formation of standing waves is possible when n is a multiple of 3, as the low-energy states about the two Dirac 
points fall then in the same sector with quantum number q — 1 regarding the Cgt, symmetry. Thus, it is possible to 
form a state confined in the nanotube by superposition of two modes with opposite longitudinal momenta and the 
same quantum number q. This is consistent with the fact that the confined states are actually found in the q = 1 
sector in the diagonalization of very large lattices of individual nanotube-graphene junctions. In the real lattice of the 
array, the confinement of the electrons in the nanotubes is only approximate, but the decay of the wave functions in 
the graphene part away from the junctions is strong enough to account for the development of the flat bands shown 
above. 



V. CONCLUSIONS 

In this paper we have studied the electronic structure of the hybrid material made of carbon nanotubes attached 
to a graphene sheet. By analyzing individual nanotube-graphene junctions, we have found the following features: 

i) Low-energy electrons in the graphene layer, with |e| ^ vp/ Ro, are scattered by the nanotube, and the probability 
of propagating into the tube is small. 

ii) High-energy electrons reaching the nanotube junction, with |e| ^ vf/Rq, are mostly transmitted into the 
nanotube. 

iii) At low energies, |e| ^ vp / Rq, and in the vicinity of the junction, r ^ Rq, there is in general a depletion of the 
density of states. 

iv) In certain nanotube geometries (armchair and (6n, 0) with n equal to a multiple of 3), there are quasi-bound 
states near e = 0, partially localized at the junction. 

We have shown that these features can be accounted for in a continuum model of the hybrid geometry. This is based 
on the Dirac fermion fields describing the electronic excitations, interacting with the curvature and the effective gauge 
field arising from the six heptagonal carbon rings at the junction. Thus, properties i), ii) and iii) are intrinsic to the 
continuum Dirac equation and universal for all nanotube geometries, while iv) depends on the relative position of the 
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six heptagonal rings and the consequent effective magnetic flux at the junction. While we have focused on the case of 
armchair and zig-zag nanotubes, it becomes clear that the continuum theory may account as well for the properties 
of junctions with other geometries. In this respect, it is likely that, by allowing for less regular distributions of the 
heptagonal rings, nanotubes with nontrivial helicity can also be attached to the graphene sheet. In a more general 
theoretical perspective, it may be interesting to analyze other discrete realizations of the 2D Dirac equation, like the 
geometry of a square lattice with one half magnetic flux per plaquette. 

We have also shown that the arrays of nanotube-graphene junctions fall into two main classes, depending on whether 
their spectra exhibit or not flat bands close to the Fermi level. The flat bands only appear in arrays made of armchair 
nanotubes or (6n, 0) nanotubes when n is a multiple of 3. On the other hand, the semiconducting behavior seems to 
be a constant in the class characterized by the presence of the flat bands, as no dispersive bands cross then the Fermi 
level. Metallic behavior of the array of junctions is possible in the other class, as shown in Fig. [5fb), though that 
behavior does not appear to be a generic trend, as illustrated by the absence of low-energy bands crossing the Fermi 
level in the other representative of the class shown in Fig. [H^d). 

In real experimental samples, it is quite likely that the arrays may be formed by junctions with nanotubes of different 
helicities. In this case, we can expect that the electronic structure of these arrays will be a mixture of the features 
already present in Figs. [5]^a)-(d). In particular, part of the electronic states will be still confined in some of the 
nanotubes, and other states will be partially localized at some of the junctions. The feasibility of using the arrays of 
nanotube-graphene junctions may depend on the possibility to tailor these hybrid structures to get specific functions. 
At this point, more input from experimental measurements on these arrays would be required, while the remarkable 
behavior predicted for these systems (localization and confinement of states, flat bands) opens good perspectives in 
the investigation of novel electronic devices. 
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